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Abstract —A novel dynamic radio-cooperation strategy is pro¬ 
posed for Cloud Radio Access Networks (C-RANs) consisting of 
multiple Remote Radio Heads (RRHs) connected to a central 
Virtual Base Station (VBS) pool. In particular, the key capa¬ 
bilities of C-RANs in computing-resource sharing and real-time 
communication among the VBSs are leveraged to design a joint 
dynamic radio clustering and cooperative beamforming scheme 
that maximizes the downlink weighted sum-rate system utility 
(WSRSU). Due to the combinatorial nature of the radio clustering 
process and the non-convexity of the cooperative beamforming 
design, the underlying optimization problem is NP-hard, and is 
extremely difficult to solve for a large network. Our approach 
aims for a suboptimal solution by transforming the original 
problem into a Mixed-Integer Second-Order Cone Program 
(MI-SO CP), which can be solved efficiently using a proposed 
iterative algorithm. Numerical simulation results show that our 
low-complexity algorithm provides close-to-optimal performance 
in terms of WSRSU while significantly outperforming conven¬ 
tional radio clustering and beamforming schemes. Additionally, 
the results also demonstrate the significant improvement in 
computing-resource utilization of C-RANs over traditional RANs 
with distributed computing resources. 

Index Terms —Cloud radio access networks; dynamic cluster¬ 
ing; joint beamforming; computing resource sharing. 


I. Introduction 

Overview: The proliferation of personal mobile-computing 
devices along with a plethora of data-intensive mobile appli¬ 
cations has resulted in a tremendous increase in demand for 
ubiquitous and high-data-rate wireless communications over 
the last few years. To cope with this challenge, the current 
trend in cellular networks is to increase the densification of 
small cells and to leverage the cooperation among multiple 
antennae and base stations (BSs). In this way, higher system 
throughput and reduced interference can be achieved via 
Coordinated Multi-Point (CoMP) transmission and reception 
techniques Q, which have been adopted in 3 GPP Long-Term 
Evolution (LTE)-Advanced. In CoMP, a set of neighboring 
cells are grouped into clusters, each consisting of connected 
BSs that share Channel State Information (CSI) and user 
signals. This scheme allows for joint processing among BSs 
that can effectively mitigate the Inter-Cell Interference (ICI) 
and thus improve the spectral efficiency. However, in current 
cellular-network architectures, physical links only exist be¬ 
tween BSs and their corresponding access network gateway 
and thus, the control signaling between BSs needed to realize 
CoMP has to travel through costly backhaul links, and often 



over a one-level higher layer in the aggregation hierarchy. 
Consequently, the latency and scarce interconnection capacity 
among BSs have resulted in limited deployments of CoMP in 
practice and, in turn, in modest BS cooperation. 

Recently, Cloud Radio Access Network (C-RAN) 0-0 
architecture has been introduced as a new paradigm for 
broadband wireless access that allows for dynamic reconfigu¬ 
ration of computing resources and provides a higher degree 
of cooperation as well as communication among the BSs. 
The fundamental characteristics of C-RAN can be summa¬ 
rized as i) centralized management of computing resources, 
ii) reconfigurability of spectrum resources, iii) collaborative 
communications, and iv) real-time cloud computing on generic 
platforms. A typical C-RAN, as shown in Fig. [U is composed 
of three main parts: 1) Remote Radio Heads (RRHs) plus 
antennae, which are located at the cell sites and are controlled 
remotely by Virtual Base Stations (VBSs) housed in a cen¬ 
tralized VBS pool, 2) the Base Band Unit (BBU) (VBS pool) 
composed of high-speed programmable processors and real¬ 
time virtualization technology to carry out digital processing 
tasks, 3) low-latency, high-bandwidth Common Public Radio 
Interface (CPRI), which connects the RRHs to the VBS pool. 

In this paper, we aim to realize the benefits offered by 
C-RANs to improve the cellular network performance via 
dynamic adaptation of radio clusters and computing resources. 
Firstly, the co-location model of the VBSs allows for their 
real-time intercommunication, thus fully enabling a coordi¬ 
nated joint transmission of RRHs that is currently practically 
constrained. In particular, control signals to realize CoMP 
between the BSs that traditionally travel via back-haul links 



can now be exchanged through the InfiniBand interconnection 
among the VBSs. A C-RAN-based radio-cooperation scheme 
would be fully dynamic and user specific, in the sense that 
we can form a virtual cluster of RRHs to coordinate their 
downlink transmissions to each of the scheduled users. In 
this strategy, each scheduled user is always the central of a 
RRH cluster, making it different from the traditional CoMP 
techniques where the RRHs are grouped into fixed and non¬ 
overlapping clusters. 

Related Works: Pioneering works on realizing the benefit 
of C-RANs have focused on the overall system architecture 
with emphasis on system issues, feasibility of virtual software 
base station stacks, performance requirements and analysis 
of optical links between RRHs and their VBSs £0, 0, 
0 . On the other hand, considerable attention has also been 
paid on cooperative communications techniques for C-RAN 
under various different objectives. For example, in 0 the 
authors propose a blind source separation strategy to mitigate 
interference in uplink C-RAN; in 0, 0 the authors consider a 
network power minimization problem. In addition, the optimal 
tradeoff design between transmit power and backhaul capacity 
is studied in (na, while the tradeoff between transmit power 
and delay performance is investigated in ifTTTl via a cross-layer 
based approach. 

In this paper, we study a dynamic radio cooperation 
technique and consider Weighted Sum-Rate System Utility 
(WSRSU) as the performance metric under a practical con¬ 
straint on computing resources at the VBS pool. Note that 
the BS cooperation for WSRSU maximization problem has 
been studied in traditional CoMP systems. However, due the 
scarce interconnection among the BSs and the lack of global 
CSI available at each BS, existing clustering and coopera¬ 
tive beamforming techniques are mostly heuristic-based (i.e., 
the clustering decision is made based on the relative signal 
strength and locations of the users, and the beamforming 
design is not adaptive to inter-cluster interference) 1X21 - 1141 . 

Our Contributions: In this paper, we propose a novel 
dynamic radio cooperation strategy for C-RANs that takes ad¬ 
vantage of real-time communication and computing-resource 
sharing among the VBSs. Unlike existing methods, our ap¬ 
proach makes the joint clustering and beamforming decision 
based on global CSI available at the VBS pool, thus being able 
to mitigate both the intra-cluster and inter-cluster interference 
in order to significantly improve the system’s performance. 
Our proposed solution dynamically groups the RRHs into 
user-specific (potentially overlapping) clusters and designs the 
downlink beamformers at each RRH in order to maximize 
the WSRSU function. In particular, within each scheduling 
interval, i.e., a time-frequency resource block, a group of 
RRHs is identified to serve each scheduled user. To realize the 
proposed Dynamic Radio Cooperation (Dynamic-RC) strategy, 
we formulate the associated optimization problem, which 
we also refer to as the Dynamic-RC problem , that aims to 
maximize the WSRSU under the transmit power constraints 
at the RRHs and the total computing-resource constraint at 
the VBS pool. Due to the combinatorial nature of the radio 


clustering process and the non-convexity of the cooperative 
beamforming design, the Dynamic-RC problem is extremely 
difficult to solve optimally in practical (polynomial) time for a 
system with a large number of users and RRHs. To overcome 
this drawback and solve the problem efficiently, our approach 
aims for a suboptimal solution with reasonable complexity. 
In particular, we exploit conic programming techniques lfT5l 
and the /1-norm reweighting approximation methods from 
Compressive Sensing which were originally proposed for 
sparse signal recovery ED, in order to quickly identify the 
optimal clustering decision and beamforming design. 

We propose an iterative algorithm to solve the Dynamic- 
RC problem. In each iteration, the clustering decision is 
temporarily fixed and a Cooperative Beamforming Design 
(CBD) problem is solved using Second-Order Cone Program¬ 
ming (SOCP) technique. The optimal beamforming solution 
obtained from the CBD problem is used to adjust the clustering 
decision via /1-norm reweighting technique. As such, the joint 
clustering and beamforming design is quickly identified and 
is adaptive to the global network condition. 

Numerical simulations are carried out extensively in various 
user distribution scenarios and demonstrate that our proposed 
low-complexity Dynamic-RC strategy significantly improves 
the WSRSU performance over conventional radio clustering 
and beamforming schemes. Furthermore, the results also show 
the great potential gains of C-RANs using our Dynamic-RC 
strategy over distributed RANs in terms of computing resource 
and transmit power utilization. 

Paper Organization: The remainder of this paper is or¬ 
ganized as follows: in Sect. [Ill we present the considered 
system model and formulate the problem under study; in 
Sect, [nil we discuss the analysis and solution to the coop¬ 
erative beamforming design problem with a fixed clustering 
decision; in Sect. uyi the Dynamic-RC strategy via dynamic 
radio clustering and beamforming design is solved via our 
proposed iterative algorithm; simulation results are illustrated 
in Sect. [V| and, finally, Sect. [VXl concludes the paper and points 
to future work. 

II. System Model and Problem Formulation 

In this section, we firstly introduce the system model 
of the considered downlink C-RAN system and discuss the 
computing-resource constraint. The proposed dynamic radio 
cooperation strategy is then formulated as a joint clustering 
and beamforming design problem. 

A. System Model 

We consider a multi-user, multi-cell C-RAN downlink 
system, where each cell has one RRH that connects to a 
common VBS pool via high-capacity backhaul links. Let 
1Z = {1, 2,..., R} be the set of RRHs and U = {1, 2,..., U} 
be the set of active users in the system. We assume that each 
RRH r has N r antennae while, realistically, all the users are 
equipped with only a single antenna. Note that the solutions 
proposed can be trivially extended to the multi-antenna-user 
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case. The RRHs cooperate with each other to form virtual user- 
specific clusters, i.e., each RRH cluster is formed for a sched¬ 
uled user, while each RRH can be part of multiple clusters. 
Hence, the number of virtual clusters is equal to the number 
of scheduled users in the system, which may be smaller than 
the number of total active users. Let S = \u G U, r G 1Z} 
be the clustering decision, in which s r u is a binary variable 
equal to 1 if RRH r is selected to serve user u, and 0 
otherwise. Consequently, let V u = {r 1Z\s r u = 1} denote 
the serving cluster of user u. We consider the system in a 
single time-frequency resource block, which is considered to 
be spatially reused across all the users. As such, each RRH can 
simultaneously serve at most N r users; otherwise, the users 
will suffer from intra-cluster interference. 

We assume that each user has a single traffic flow that is 
independent of all other users’ flows. Baseband signals for user 
u and the corresponding downlink beamforming information 
after being processed at the VBS pool will be transported to 
all the RRHs in the serving cluster V u . In each scheduling 
slot, all the RRHs in V u will jointly transmit the normalized 
symbol x u G C of unit power to user It is assumed that 
the signals for different users are independent from each other 
and from the receiver noise. Now, let G qn t xi 
linear downlink beamforming vector at RRH r corresponding 
to user u and W = {w^ \\/u £U,r £1Z} denote the network 
beamforming design. Note that W also implies the scheduling 
decision, i.e., user u is not scheduled for the current time- 
frequency slot if = 0,Vr G 1Z. In the current scheduling 
slot, the received signal y u G C at user u is, 

y u = h u™u x u + h u W u' X u' + z u, ( 1 ) 

r£V u u',u'^£u r'£V U / 

v ---' '-v-- 

desired signal interference 


infrastructure of a datacenter. Each VBS performs baseband 
processing for a certain set of users, and by leveraging virtual¬ 
ization technology, these VBSs can flexibly share the common 
computing resource of the physical server pool. Recently, 
the implementation of software VBSs on General-Purpose 
Platform (GPP) has been realized (see, for example, 0, 
0). Profiling results on these systems have revealed that 
the utilized computing resource at a VBS is an increasing 
function of the accumulated data rates processed by that VBS. 
Therefore, it is reasonable to argue that the total computing- 
resource capacity of the VBS pool places a cap on the 
total data rates of the users in the network. In general, the 
computing-resource capacity of the VBS pool can be modeled 
as a multi-dimentional vector representing the capacities of the 
CPUs, memory, and network interfaces. However, for the ease 
of analysis, we only consider scalar computing capacity in this 
paper. In particular, let C denote the total computing capacity 
in the VBS pool that can be flexibly shared among all the 
VBSs. The computing-resource constraint on the accumulated 
data rate of all the users in the system can be expressed as 



where R u is the data rate of user u given in (0 and T(.) is 
an increasing function specifying the relationship between the 
utilized computing resource and the accumulated user data 
rat^E It should be noted that for a traditional system with 
distributed computing resource at the RRHs, the accumulated 
data rate processed at each RRH r will be subject to the per- 
RRH computing-resource constraint C r < C , i.e., 

<c r ,vren. (5) 


where G C lxNr is the channel coefficient vector from 
RRH r to user u, z u is the zero-mean circularly symmetric 
Gaussian noise denoted as C./V(0,cr 2 ). For simplicity, let 

= E h u W u' and = Vu,u = E h n w n- With 

r'EV u / rEV u 

this position, the received Signal-to-Interference-plus-Noise 
Ratio (SINR) at user u is, 


Ifu 


l^l 2 


( 2 ) 


B. Joint Clustering and Beamforming Problem Formulation 

Our objective is to maximize the WSRSU under the transmit 
power constraint at each RRH and the total computing- 
resource constraint at the VBS pool. It is assumed that the 
capacity of the front-haul links connecting RRHs to the VBS 
pool is sufficiently provisioned to accommodate peak-capacity 
demand. Our proposed dynamic radio cooperation strategy 
involves finding the optimal clustering decision <S* and the 
optimal beamforming design W*, and can be formulated as, 


LL i CL -f— LL 

Thus, under the clustering decision S and the beamforming 
design W, the Shannon transmission rate of user u can 

(S*,W*)= argmax 

reiz,ueu 

V ( i» R u (s, W) 

uEU 

(6a) 

be calculated as R u (S, W) = yBlog 2 (1 + u ), in which 

B [Hz] is the channel bandwidth and [0,1] account for 

s.t. 

Il w «ll2 < p r,Vr G K, 

uEU 

(6b) 

the spectral and the coding efficiencies, respectively. Unless 
otherwise stated, for notation simplicity in the subsequent 


ll w «H 2 < s u P r, 

(6c) 

analysis we will assume B = r\ = p = 1 and consider the 
normalized rate (bits/s/Hz). Hence, the rate R u simplifies to, 

r> / c i i \ 


Y R u(s, w) <n, 

uEU 

V s « ^ Nr, S T U G {0, 1} , 

(6d) 

(6e) 


(3) 


uEU 


Computing resource constraint The VBS pool consists of i The realization of r( . } can be obtained by carefully profiling the VB Ss 
a set of interconnected VBSs hosted in the physical-server at different level of offered load in a practical C-RAN implementation. 
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where q u ,u eU, is the utility marginal function corresponding 
to user u, which can represent the user-specific Quality of 
Service (QoS) or priority in the system, P r [W] is the per- 
RRH transmission power constraint and Cl = argr(C). 
Constraint (l6cl) indicates the coupling between the assignment 
variable s r u and the beamforming vector w£, i.e., wjj = 0 
when s r u = 0. We refer to © as the dynamic radio cooperation 
(. Dynamic-RC ) problem. In fact, this is a Mixed-Integer Non- 
Linear Program (MINLP), which is intractable in practical 
time. Specifically, even when the binary variables s r u are fixed, 
solving for is still NP-hard. 

Given a large number of variables that scales linearly with 
the number of users and RRHs in the system, finding a low- 
complexity, suboptimal solution is highly desirable. To this 
end, we firstly solve the Cooperative Beamforming Design 
(CBD) problem with given clustering decision S , and propose 
a low-complexity iterative algorithm to solve the Dynamic- 
RC problem to a local optimum. Specifically, in Sect. [Till 
we will transform the CBD problem into a SOCP with a 
fixed clustering decision, and will take advantage of the 
existing efficient SOCP algorithms. The Dynamic-RC problem 
will then be solved in Sect. U3 using the iterative /1-norm 
reweighting technique, which solves the CBD problem and 
updates the clustering decision in each iteration. 


III. Cooperative Beamforming with Fixed 
Clustering Decision 


In this section, we consider the problem of Cooperative 
Beamforming Design (CBD) for a given radio clustering deci¬ 
sion S. In particular, for given {s^} satisfying constraints (l6el) . 
we need to find the optimal downlink beamformers {w^} by 
solving the CBD problem below, 


max V' q u R u ( S , W) 

wr,rG7 Z,uEU J 

(7a) 

ueU 


S-t. J2 Il w «ll2 < ^r,Vr G K, 

(7b) 

uEU 


y^Ruis, w) < n. 

(7c) 

u£U 



Observe that the rate functions R u ’s appear in both the 
constraint and objective of ©, making the problem diffi¬ 
cult to deal with. To decouple this problem with respect to 
(w.r.t.) Ru’ s, we remove the constraint (TTcl) and consider the 
relaxed-CBD problem with constraint (1751) only. The solution 
{w^} of the relaxed- CBD problem will be verified against 
constraint G3 so to finally obtain the solution of the original 
CBD problem by solving an additional feasibility problem. In 
the following subsections, the relaxed-CBD first and then the 
feasibility problem will be addressed sequentially. 


A. Relaxed-CBD Problem 

The relaxed- CBD problem is rewritten from © without the 
computing-capacity constraint G3, and is cast as follows, 

max 'V' q u R u (S, W) (8a) 

wi;,rG7^,uG^/ ^ J 
ueU 

s.t. Il w «ll 2 < JV,Vr € (8b) 

udlA 

This is in fact a weighed sum-rate maximization problem, 
which is widely known to be NP-hard. Our approach aims for 
a local solution using a low-complexity algorithm designed 
by effectively exploiting the techniques of SOCP0 In order 
to use the efficient algorithms developed for SOCP, one 
must reformulate the problem into the standard form that the 
algorithms (e.g., those proposed in 03) are capable of dealing 
with. Firstly, from ©, objective function (l8at is rewritten as, 

Y, (s, w) = yy i og2 (i+ 7 «) 9 ". (9) 

Now, by introducing the variables £ u ’s, u eU, we can recast 
the relaxed-CBD problem in ^ as, 


max TT t u (10a) 

wT ,rElZ,u£U 

uEU 

s.t. 7 „ > -lyueU, (10b) 

£lKll 2 ^A,Weft, (10c) 

uEU 

which stems from the fact that constraints (II Obi) are active at 
the optimum. We now have the following Lemma. 

Lemma 1. Let where f r u is the phase rotation 

such that the imaginary part of equals to zero, Mu G 

U,r G IZ. If is optimal to m, then is also optimal 

Proof: We can represent as = |h^wjj| e j6 ^. 

By choosing f r u = -0 r u , we have = h^w^e- 7 ^ = 

|h^wjj|. Recall y u given in ©, it is straightforward to verify 
that substituting by w£, Mu E U,r E 1Z, into (ITOl) will 
result in the same objective function and constraints. Thus, if 
is optimal then is also optimal. ■ 

Using Lemma 1, we can restrict ourselves to the beamform¬ 
ers in which > 0, Mu G U, r G V n , where each product 

has a non-negative real part and a zero imaginary part. Notice 
that constraint (II Obi ) is equivalent to 

- ——2 - >ti /,u -i,VueZZ (li) 

E \*u,u ’\ 2 + V 2 ~ 

u' ClIA,u'^u 

which can be recast as, 




and y2 1^“-“' 


( 12 ) 

2 +a 2 <p u ,Wu€U, (13) 


2 Second-Order Cone Problems (SOCP) are convex-optimization problems 
in which a linear function is minimized over the intersection of an affine set 
and the product of second-order (quadratic) cones. 
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by introducing the slack variables /3 u ’s and due to the fact 
that both constraints (IT2l) and (fl3l) are active at the optimum 
of problem (fTOl) . It can be verified that (IIPel) and (fl3l) 
follow the Linear Programming (LP) constraint expression 
with generalized equalities/inequalities, which can be directly 
written as Second-Order Constraints (SOCsjl |33). To deal 
with the non-convex constraint m, we further exploit the 
sequential parametric convex-approximation approach in fHOl 
to approximate (IT2l) as convex as presented in the following. 

Firstly, (Eli can be rewritten as, 

(14) 

& + 1 >ti /a “,VueW. (15) 

Observe that, for a given <p u , we have 

( 16 ) 

which follows the inequality of arithmetic and geometric 
means of 4> u Pu an d The equality in (fl6l) is achieved 

when (j) u = and we get the equivalent form of 

constraint Eli as, 

(17) 

20 u 2 

Furthermore, without loss of generality, we scale q u ’s in (TTab 
such that q u > l,\/u G U to make tl/ qu become concave. 
Thanks to the concavity of t u ’s, we can adopt the results 
in ED to replace the right side of (fl 51 ) by its iterative first- 
order approximation as, 

+ (18 ) 

Qu 7 ' 

where t y u J denotes the value of t u in the previous iteration. 
From (fl3l) . (fT7t . and (fl8l) . the re/a*ed-CBD optimization 
problem in ([8]) can be finally recast as, 


max TT t u 

WT ,r£7Z,u€:U 

ueu 

(19a) 

11 11 2 < PrXr e n, 

(19b) 

uCQA 


dmi,C3.C9. 

(19c) 


Notice that the objective function and all the constraints 
in © admit SOC representation (see CD, 03). Conse¬ 
quently, the resulting problem in (IT9b is a SOCP, which can 
be solved efficiently and very fast using standard solvers such 
as CPLEX (H or MOSEK l20l . 

B. CBD Feasibility Problem 

Here, the solution of the relaxed-CBD problem ([8]) which 
was obtained via solving the equivalent SOCP problem in (fl9l) . 
will be verified against the computing-capacity constraint 
in (I7c1) to obtain finally the beamforming solution of the 
original CBD problem cast in 0). 

3 In a SOC representation, the hyperbolic constraint ab > c 2 , with a, b > 0, 
is equivalent to 11 [(a — b) 2c] T | |2 < a + b. 


Suppose that W is the beamforming solution of prob¬ 
lems ([8]). If W satisfies the computing-resource constraint 
G3, i-e., E Ru [S, W) < ft, then W is also the optimal 

ueu ' ' 

solution of 0). In this case, the WSRSU is limited by 
the per-RRH power budget only, and not by the computing- 
resource capacity of the VBS pool. On the other hand, when 
the computing-resource constraint is violated, we need to 
selectively drop the rates of some users. This can be done 
via a greedy algorithm that keeps dropping the users that 
have the smallest marginal utility function q u from the current 
scheduling interval until the total data rate of all the scheduled 
users satisfies the computing-resource constraint. Since the 
optimal bearmformer design W is jointly calculated for all 
users, dropping the rates of some users requires recalculating 
the beamformers of all the RRHs. 

Let { R * > 0 ,u eU} be the user rates obtained after the 
greedy-user-rate-dropping process is applied; the beamformer 
design W that achieves these rates can be obtained via solving 
the feasibility problem given below, 


find 

s.t. 


{ W u}, M 

eu,r eV u 

(20a) 

X H w « 

|| 2 2 < P r yrell, 

(20b) 




E 

lip I 2 

' 1 , 2 —- > 7 ;,vm€w, 

(20c) 




u' CQA ,u' 

where 7 * = 2 R « — 1 . 

The feasibility problem in (f20l) is not convex', however, by 
exploiting its special structure, we can transform this problem 
into a SOCP form, which can be solved efficiently. The trans¬ 
formation is presented as follows. Firstly, let w r be the long 

column vector such that w r = (w^) T , (w 2 ) T , ...(w[}) T , 

Vr G 7 Z. Constraint (I20bl) can be rewritten in a SOC form as 

l|w'| a <V^,Vrefc. ( 21 ) 


Furthermore, (I20cl) is equivalent to 


1 

7 u 


|^«| 2 > X | 2 +* 2 .Vr eU. (22) 


u'EU 


Since > 0, as we considered previously, we can take 

the square root of both sides in (l22l) . which yields, 






(23) 


It can be seen that (f23l) follows the SOC form; hence, 
using (fTTI) and (l23l) , we are now ready to recast the feasibility 
problem in (l20l) in the standard SOCP form as follows, 


find (wy,«6W,rg V„ 

(24a) 

iiw i 'n 2 < yy r ,\/r e n, 

(24b) 


(24c) 

V 7 « 
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The solution W* for (t24l) can be obtained using standard 
SOCP techniques such as the interior-point methods EH or 
the SOCP solvers (e.g., CPLEX, MOSEK). In summary, the 
optimal beamformer design of the CBD problem in © for 
a given radio clustering decision S can be obtained by the 
procedures described in Algorithm [l] 


Algorithm 1 Cooperative Beamformer Design (CBD). 

(1) Solve the SOCP problem in (IT9l) to find W 

(2) Verify constraint (fTcb 

• If T, Ru (S, W) < return W*=W. 

ueu ' I 

• Otherwise: Drop users’ rates using the greedy algorithm 

- Repeat : Update R u > = R u ' (jS, — r, 

where r is small decreasing step and 

q u '=min{q u : q u > 0,u G U}. Go to the next 
user when R u f ^5, = 0. 

- Until : Ru [5,w] < Q 

ueu v J 

- Solve feasibility problem (1201) , get W*. Return 


IV. Joint Dynamic Radio Clustering and 
Beamforming Design 


In the previous section, the CBD problem has been trans¬ 
formed into an equivalent SOCP form. As a result, our con¬ 
sidered Dynamic-RC problem in © can also be transformed 
into a Mixed-Integer SOCP (MI-SOCP) problem with binary 
variables s^’s. In a network with U users and R RRHs, there 
are 2 UR possible clustering patterns. The optimal solution to 
the clustering decision can be found via exhaustive search 
or using standard global optimization solvers. However, these 
approaches usually have a complexity growing exponentially 
with the problem size, which is not a practical approach. 
Hence, in this section, we present a method to solve the 
Dynamic-RC problem given in © by iteratively solving the 
CBD problem using Algorithm [T] In particular, we take 
advantage of the /1-norm reweighting technique to adjust the 
approximation of the clustering variables after each iteration. 

Firstly, given the relationship of s r u and w£, we can repre¬ 
sent s r u by /0-norm expression of w r u as follows, 


s 


r 

u 



(25) 


The above expression allows us to leverage the /1-norm 
reweighting technique, which has been effectively applied in 
the literature to approximate the /0-norm lfl6lK i.e., ||x || 0 « 
Y^PkXk, where x C M n and pi, p 2 ,p n are positive weights. 

k 

With n = 1, by choosing Xk = || || 2 , we get s r u « 

Pu || W u || 2’ i n which the weight p r u is adjusted iteratively as 

Pu= „ yueu,ren, (26) 

Mla + e 

with ||H 2 obtained from the previous iteration. In d26b . the 
parameter e is a very small positive number introduced to 


provide stability and to ensure that in case || ||^ = 0, it does 

not strictly prohibit a non-zero estimate in the next iteration. 

The Dynamic-RC problem in © - given now s - can 
be rewritten as, 


max Tq u R u (5,W) 

wj;,rG7^,uGW ^ J 
uEU 

(27a) 

s.t. X HwUla < PrXr € n, 

uCiU 

(27b) 

(s, w) <n, 

uEU 

(27c) 

YsPu H W «ll 2 ^ N r- 

uEU 

(27d) 


Note that constraint (I27db can be written in SOC form as, 




< \/~N~ r , Vr G TZ. 


(28) 


Thus, the problem in (1271) is similar to the CBD problem 
in © with the additional SOC constraint (l28l) . which can be 
solved efficiently using Algorithm [l] To clarify the idea, we 
present the iterative method to solve the Dynamic-RC problem 
in Algorithm [2] below. 


Algorithm 2 Dynamic Radio Cooperation via Iterative SOCP 

(1) Initialization: set p r u = 0, Vu e U, r e 1Z 

(2) Iteration: 

a) Solve problem (l27l) with the current value of p r u using 
Algorithm [I] In particular, Step (1) in Algorithm [I] will 
solve problem (fl9l) with the additional constraint (l28t . 

b) Update the weights p\f s using the solution w^’s ob¬ 
tained in the previous step as, 

P r u = (\\™u\\l + e ) ,VueU,reK. (29) 

(3) Check convergence: Repeat Step (2) until convergence 
or the max number of iterations is reached. 


Note that RRH r is included in the serving cluster of user 
u, i.e., r G V u , if the beamformer from RRH r to user 
u , w£, is nonzero. Since p r u = 0 in the first iteration in 
Algorithm 2, the constraint (l28l) is automatically satisfied. 
Thus, initially each RRH can be selected into more than N r 
clusters. After that, the weights {p r u } are updated inversely 
proportional to the beamforming power as in (|29b . Therefore, 
among the beamformers from all the RRHs to a target user, 
those with highest powers are most likely to be identified 
as nonzero in the next iteration. This allows for successive 
better estimation of the clustering decision, i.e., identifying 
the nonzero beamformers from RRHs to users. As will be 
shown later in our simulation results, the beamforming powers 
quickly converge within a few iterations. 

Complexity analysis: The computational complexity of Al- 
gorithm [2] mainly lies in Step (2a) where a SOCP problem 
is solved. Assuming the same number of antennae N r on the 
RRHs, the total number of variables in this SOCP problem 
is URN r , where U and R are the numbers of users and 
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RRHs. Thus, the computational complexity of the interior- 
point method to solve such a SOCP problem is approximately 
O (jJJRN r ) 3 ‘ 5 ^j lETI . This is significantly advantageous for 
a large network compared to the optimal design using ex¬ 
isting solvers, which are characterized by a prohibitively 
exponential-time complexity. 

Furthermore, in practical networks, a RRH r should not be 
included in the serving cluster of user u if r is very far away 
from u. Assuming a network of hexagonal cells, we can pre¬ 
select only the 7 RRHs having strongest channel coefficients to 
user u to be the candidate serving cluster of user u, denoted as 
C u . After the pre-selection process, Algorithm 2 will identify 
the optimal serving cluster V u within the subsets of C u . This 
can significantly reduce the complexity of Algorithm 2 to 
O 7UN r ) 3 ' 5 ^. We adopt the pre-selection of serving cluster 
candidates in the simulation and numerical results show that 
this approach performs very close to the optimal solution. 

V. Performance Evaluation 

In this section, simulation results are presented to evaluate 
the performance of our proposed Dynamic-RC algorithm. We 
consider a network of hexagonal cells with a RRH in the 
center of each cell. The neighboring RRHs are separated 
1 Km apart from each other. We assume that all the wireless 
channels in the system experience block fading such that 
the channel coefficients stay constant during each scheduling 
interval but can vary from interval to interval, i.e., the channel 
coherence time is not shorter than the scheduling interval. We 
assume that all the RRHs have the same number of transmit 
antennae N r and transmit power budget P r . The channel 
coefficients are calculated following the path-loss model, given 
as L [dB] = 148.1 + 37.6 log 10 d[ km ], and the log-normal 
shadowing variance set to 8 dB. In addition, it is assumed 
that the channel bandwidth B is 10 MHz, is reused across all 
the users, and the noise spectral density is —100 dBm/Hz. 

WSRSU performance: Firstly, we consider a system with¬ 
out the computing-resource constraint and evaluate the perfor¬ 
mance of the four radio cooperation algorithms below. 

• Optimal: The WSRSU of the optimal scheme is obtained 
by using the solver MOSEK to solve the equivalent MI- 
SOCP presentation of problem ©. 

• Dynamic-RC : Our proposed dynamic radio cooperation, 
where the solutions are obtained from our iterative, low- 
complexity Algorithm [2] 

• CVSINR : A downlink cooperation scheme proposed in 
m where the cluster for each user is formed heuris- 
tically based on the relative signal strength and the 
clustered virtual SINR (CVSINR) algorithm is used to 
design the beamforming vectors. 

• Greedy : A greedy clustering algorithm proposed in lfl4ll . 
which solves an equivalent set covering problem to 
select the set of non-overlapping base station clusters. 
This scheme uses zero-forcing as the criterion to design 
beamformers and a greedy algorithm is used for user 
scheduling. 



Scenario 1 (Uniform) Scenario 2 (Uneven) Scenario 3 (Extremely Uneven) 

Light Medium |||j Heavy 

Fig. 3. Different user distribution scenarios: Scenario 1 (uniform) with all 
medium (loaded) cells; Scenario 2 (uneven), light and heavy (loaded) cells are 
intermixed together; Scenario 3 (extremely uneven), heavy cells are grouped 
together, and the heavy cell group is surrounded by light cells. 


We evaluate the four schemes above in a network of 16 
cells with three different user distribution scenarios as shown 
in Fig. 0 In particular, Scenario 1 consists of all medium 
(loaded) cells where users are distributed uniformly over all 
the cells; Scenarios 2 and Scenarios 3 consist of light and 
heavy (loaded) cells, however the heavy cells are intermixed 
with light cells in Scenarios 2 to represent micro-tidal effect 
while they are grouped together in Scenarios 3 to represent the 
macro-tidal effect. In our simulation, we perform 500 drops, 
in each drop 32 users are placed randomly in the network with 
1 user in a light cell, 2 users in a medium cell and 3 users 
in a heavy cell. The utility marginal functions q u ’s are chosen 
randomly such that 0 < q u < 1, \/u G U. 

Fig. [2]-(a), (b), (c) plot the WSRSU performance of the 
four considered radio cooperation schemes in Scenario- 1, 2, 
3, respectively. It can be seen that our proposed Dynamic- 
RC scheme and the Optimal scheme significantly outperform 
the CVSINR and Greedy schemes in all three scenarios. This 
is because the heuristic clustering of the RRHs in the later 
two schemes is suboptimal, plus their beamforming design 
algorithms only aim to minimize the intra-cluster interference 
but not the inter-cluster interference. On the other hand, our 
proposed Dynamic-RC scheme takes into account the global 
network condition that is available at the VBS pool, which 
provides better clustering decision and beamforming design. 
Compared to the optimal scheme, our proposed Dynamic- 
RC strategy via Algorithm [2] shows a small loss in WSRSU 
performance but has a significant advantage in reducing the 
execution time. In fact, in our simulation for the considered 
system configuration (U= 32, R= 16), MOSEK solver takes 
more than 100 s to obtain the optimal solution of the MI- 
SOCP problem, while each iteration in Algorithm [2] takes less 
than a second and the algorithm overall converges within 15 
iterations. 

Impact of Maximum Cluster Size: Fig. S-(a), (b) plot the 
CDF of average user rate (w.r.t. 32 users) achieved by 
Dynamic-RC scheme with different choices of the maximum 
cluster size, Vmax • In each case, only V m ax RRHs having 
the strongest channel coefficients to a user are chosen to 
be the candidates of that user’s serving cluster. This pre¬ 
selection is done before running Algorithm 2 to finally find 
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500 





(a) (b) (c) 

Fig. 2. Weighted Sum-Rate System Utility (WSRSU) of a C-RAN downlink system using different radio cooperation schemes, evaluating on three different 
user distribution scenarios. (a)-Scenario 1, (b)-Scenario 2, (c)-Scenario 3. 




Fig. 4. CDF of Average User Rate obtained by Dynamic-RC scheme with 
different numbers of the maximum cluster size. (Pr = 10 dBm) 


the best serving cluster for each user. When Vmax = 1, there 
is no cooperation among the RRHs. The results in Fig. 0 
(a), (b) are obtained by performing 500 drops on Scenario 
2 and Scenario 3, respectively, with Pr =10dBm. The utility 
marginal functions are updated in each drop according to the 
proportional fairness criterion, i.e., q u = 1 /R u where R u 
is the long term average data rate for user u E U. It can 
be seen that the improvement in average user rate due to 
larger cluster size in Scenario 3 (macro-tidal effect) is greater 
than that of Scenario 2 (micro-tidal effect). For example, the 
dynamic cooperation scheme with V ma x = 3, 5, 7 provides 
130%, 137%, and 138.6% gain, respectively, for the 60th- 
percentile average user rate over the non-cooperation scheme 
{Vmax = 1), in Scenario 2; while the corresponding gains in 
Scenario 3 are 145%, 159%, and 162%, respectively. Although 
not included here due to space limitation, we observe that 
when Vmax exceeds 7 cells, the additional gain is negligible. 

Benefits of Computing Resource Sharing: To evaluate the 
impact of the computing-resource constraint on the system 
performance, Fig. [5] compares the WSRSU performance of our 
considered system with the centralized computing-resource 


constraint, as expressed in (0]), versus a conventional system 
with a distributed computing-resource constraint, as expressed 
in (0. In particular, we consider a network of 4 cells with 2 
users in each cells in random locations and q u ’s are chosen 
randomly. For a fair comparison, we set arg T(C) to 400 Mbps 
and argr(C r ) to 100 Mbps, and ran the Dynamic-RC scheme 
in Algorithm [2| on both systems. Note that, in this setting, 
each of the 4 RRHs in the distributed system is provisioned to 
process maximum 100 Mbps of user baseband traffic at a time, 
while in the centralized system the VBS pool is provisioned 
to process maximum 400 Mbps baseband traffic at a time. We 
say that the computing resource is saturated in each system 
when the achieved sum-rate (SR) of all the users reaches the 
maximum provisioned processing traffic rate. As the transmit 
power increases, observe in Fig. 0 that the computing capacity 
of the VBS pool in the centralized system saturates earlier 
than the total computing capacity of the distributed system 
does (when the computing capacity is saturated at all the 
RRHs). In fact, the WSRSU and SR of the distributed system 
saturate almost at the same time while the WSRSU of the 
centralized system continues to increase after the saturation 
point (of the SR), and is significantly higher (up to 250% 
gain) than that of the distributed system. This demonstrates 
the great potential gains of C-RANs using our Dynamic-RC 
scheme over the conventional distributed RANs in terms of 
WSRSU, computing resource and transmit power utilization. 

Convergence Behavior of Algorithm 2: Fig. [6] illustrates 
the convergence behavior of Algorithm [2] in identifying the 
RRH cluster for a user. We choose randomly a user u * and 
monitor the beamforming powers from the 7 candidate RRHs 
for it’s serving cluster. The evolution of the beamforming 
powers in dBm/Hz from these RRHs to user u* 9 calculated 
as || || 2, r = 1,...,7, is shown in Fig. 0 Observe that 

after the 4-th iteration, only the beamformers from RRH 1 
and RRH 4 maintain a non-trivial power, while the rest are 
forced to almost zero. In this case, the optimal serving cluster 
of user u* is identified to be V u * = {RRH 1, RRH 4} within 
only a few iterations, which demonstrates the efficiency of 
our proposed Dynamic-RC algorithm in quickly making the 
clustering decision and beamforming design. 
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Fig. 5. Downlink WSRSU of a C-RAN system with centralized computing 
resource ( Centralized ) versus that of a traditional system with distributed 
computing resources (Distributed). 



Fig. 6. Convergence behavior of the beamforming power at 7 candidate 
RRHs for the serving cluster of an example user. 

VI. Conclusions and Future Work 

We proposed a novel dynamic radio cooperation strategy 
for Cloud Radio Access Networks (C-RANs) that takes ad¬ 
vantage of real-time communication and computing-resource 
sharing among Virtual Base Stations (VBSs). The underlying 
optimization problem was formulated as a mixed-integer non¬ 
linear program, which is NP-hard. Our approach transforms 
the original problem into a Mixed-Integer Second-Order Cone 
Program (MI-SOCP) that is efficiently solved using a novel 
low-complexity, iterative algorithm. Simulation results showed 
that our low-complexity algorithm provides close-to-optimal 
performance in terms of weighted sum-rate system utility 
while significantly outperforming conventional radio cluster¬ 
ing and beamforming schemes. 

Future Work: The goal of our future work is to address 
the system-related issues and evaluate the feasibility and 
performance of the proposed strategy in a practical system. 
In fact, we are implementing a C-RAN testbed which consists 
of an open-source LTE platform OpenAirlnterface running on 
a general-purpose desktop server to realize the VBS pool, and 
a number of USRP B210/X310 boards to realize the RRHs. 
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